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ABSTRACT 

Nucleic acids are useful for biomedical targeting and sensing applications in which the molecular 
environment is different from that of a dilute aqueous solution. In this study, the influence of var- 
ious types of mixed solutions of water and water-soluble organic compounds on RNA was investi- 
gated by measuring the catalytic activity of the hammerhead ribozyme and the thermodynamic 
stability of an oligonucleotide duplex. The compounds with a net neutral charge, such as poly( eth- 
ylene glycol), small primary alcohols, amide compounds, and aprotic solvent molecules, added at 
high concentrations changed the ribozyme-catalyzed RNA cleavage rate, with the magnitude of 
the effect dependent on the NaCI concentration. These compounds also changed the thermodynamic 
stability of RNA base pairs of an oligonucleotide duplex and its dependence on the NaCl concentra- 
tion. Specific interactions with RNA molecules and reduced water activity could account for the 
inhibiting effects on the ribozyme catalysis and destabilizing effects on the duplex stability. The salt 
concentration dependence data correlated with the dielectric constant, but not with water activity, 
viscosity, and the size of organic compounds. This observation suggests the significance of the 
dielectric constant effects on the RNA reactions under molecular crowding conditions created by 
organic compounds. 

© 2014 The Authors. Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies. This 
is an open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Cellular activities occur in the presence of large amounts of 
background proteins, nucleic acids, metabolites, osmolytes, and 
cellular organelles. There are many reports suggesting that the 
molecular environment created by high concentrations of back- 
ground composition is an important driving force for several cellu- 
lar activities, including cell growth, metabolism, tumor generation, 
and age-related disease susceptibility [1-3]. The thermodynamics 
and kinetics of the chemical reactions occurring in a living cell 
are influenced by coexisting molecules through direct and indirect 
interactions. However, the structural stability of nucleic acids in a 
relevant context of living cells remains elusive. For the quantitative 
analysis of the influence of such non-dilute aqueous conditions, the 
mixed solutions of water and water-soluble organic compounds 
are used [4.5]. Neutral compounds are often employed to mimic 
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the intracellular conditions where macromolecules are concen- 
trated and the solution properties are different from pure water 
(e.g., reduction in the water activity). The non-dilute condition is 
created not only in the medium of living cells but also on the sur- 
face of a biosensor ]6-9]. Poly( ethylene glycol) (PEG) is one of the 
most commonly used cosolutes for creating a molecularly crowded 
condition. It is reported that the molecular crowding by PEG 
increases the activity of several protein enzymes, attributed to an 
enhanced folding stability and increased binding with substrates 
caused by the excluded volume effect ]10]. The polar and charged 
nature of nucleotides makes the molecular crowding effects on 
nucleic acids more complicated. For example, the thermal stability 
of oligonucleotide duplexes is not significantly changed in solu- 
tions containing a large PEG molecule, but decreases in solutions 
containing a small PEG molecule due to specific interaction with 
single strands and reduction in the water activity ]n-14]. 

RNA secondary structures can be predicted by using the ther- 
modynamic equilibrium parameters based on the nearest-neigh- 
bor model [15[. Since the structural stability of nucleic acids is 
strongly dependent on the salt concentration, empirical expres- 
sions have been developed for evaluating the salt concentration 
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dependence of the thermodynamic stability of oligonucleotide 
structures [16,17]. Intriguingly, the salt concentration dependence 
of the stability of DNA duplexes changes in mixed solutions with 
PEG [11,18-20]. A detailed study of the nature of the salt concen- 
tration dependence is important for understanding the factors 
responsible for the molecular crowding effects and for predicting 
nucleic acid structures under non-dilute conditions. In this study, 
we investigated the hammerhead ribozyme activity and the ther- 
modynamic stability of an RNA oligonucleotide duplex under the 
non-dilute conditions using mixed solutions of water and water- 
soluble neutral compounds. The hammerhead ribozyme is one of 
the smallest RNA catalysts that utilize metal ions to hydrolyze a 
target RNA sequence. Since divalent ions promote a rapid cleavage, 
many studies have been performed using Mg^* [21 ]. The Mg^* ions 
contribute to tertiary folding, and the binding is site-specific in 
some cases ]22,23]. It is a remarkable feature that the ribozyme 
motif exhibits the catalytic activity when the reaction buffer con- 
tains monovalent metal ions at high concentrations instead of 
Mg^* ]24,25]. Monovalent ions do not supposedly play an essential 
role in the catalytic step, but predominantly assist in the folding of 
an inactive conformation with a disordered catalytic core into the 
catalytically active conformation [26]. A molecular dynamics sim- 
ulation study has demonstrated that the binding patterns of Na* 
ions necessaiy for the formation of a catalytically active conforma- 
tion are not as rigid as those of Mg^* [27]. The relatively nonspe- 
cific electrostatic interactions with Na* ions are useful for 
assessing the influence of non-dilute conditions on RNA folding 
and its salt concentration dependence. The present study measures 
the rates of the hammerhead ribozyme-catalyzed RNA cleavage 
mediated by NaCl and the thermodynamic stability of a short 
RNA duplex, using aqueous solutions containing a high concentra- 
tion of neutral additives (referred to as cosolutes), such as PEG, 
small primary alcohols, amide compounds, and aprotic solvent 
molecules. This study illustrates correlations of the dielectric con- 
stant with the magnitude of the salt concentration dependence of 
the cleavage rate and the base-pair stability. The results are useful 
when constructing a prediction model for nucleic acid structures 
under non-dilute conditions and also when modulating the struc- 
tural stability by changing the solution composition. 

2. Results 

2.1. The hammerhead ribozyme reaction mediated by NaCl 

RNA sequences used for the ribozyme-catalyzed RNA cleavage 
are the minimum sequence motif of the trans-acting hammerhead 
ribozyme and the substrate strand labeled with a fluorophore, as 
shown in Fig. lA. The substrate RNA is site-specifically cleaved at 
a phosphodiester bond located in the catalytic site at the junction 
of three stems forming the pairwise coaxial stacking 123,26]. We 
confirmed that the reaction rapidly proceeded in solutions contain- 
ing Mg^* at several to several tens of millimolar, with a half-life of 
less than 1 min; but there was no reaction without adding the 
divalent metal ion to the solution. The cleavage reaction also pro- 
ceeded in solutions containing NaCl at high concentrations instead 
of MgCl2. Although the reaction rate with NaCl was much slower 
than that with MgCb, the site-specific substrate cleavage was 
observed (Fig. IB). It was previously proposed that the ribozyme 
structure in the absence of divalent metal ions is extended with 
a disordered core and no pairwise coaxial stacking [26[. The bind- 
ing of Na* ions converts to the more condensed, catalytic active 
form in equilibrium with inactive species. 

The observed reaction rate constant, k, was greatly changed 
depending on the concentration of NaCl, e.g., 2.1 h"' with 1.5 M 
NaCl, 0.31 h-' with 1 M NaCl, and 0.01 1 h^^ with 0.5 M NaCl. A lin- 
ear correlation was obtained between the common logarithm of 



the rate constant (log k) and that of NaCl concentration (log [NaCl]) 
(Fig. SIA, in the Supporting information), suggesting that the rate- 
limiting step remains unchanged within the concentration range. 

2.2. RNA cleavage rates in the presence of PEG 

The reaction with NaCl was investigated using the solutions 
containing PEG at 20 wt%, comparable to the amount of biomole- 
cules present in cells 128,29]. The effects of PEG with different 
average molecular weights (gmol^^) ranging from 2x10^ to 
2 X 10^ (PEG200, PEG600, PEG2000, PEG8000, and PEG20000) 
were compared. As observed in a previous study on the reaction 
with Mg^* 130], the RNA cleavage rate was increased by the addi- 
tion of PEGs of different molecular weights (Fig. IB). The rate accel- 
eration was greater with a larger PEG, but saturated at a molecular 
weight exceeding 2 x 10^ (Fig. IC). Remarkably, the PEG-induced 
rate acceleration was more pronounced with moderate rather than 
high concentrations of NaCl, e.g., 35-fold acceleration of the rate 
constant with 0.5 M NaCl compared to a maximum acceleration 
of 1.7-fold with 1.5 M NaCl. 

Linear dependences between log k and log INaCl] with a 
correlation coefficient greater than 0.995 were obtained, as 
shown in Fig. ID. The slope of the plot was changed in the 
PEG-containing solutions: the slope was 4.9 without cosolutes, 
and the slope became less steep as the size of the PEG increased, 
e.g., 2.1 with PEG8000. The results indicate that the PEG-induced 
rate accelerations, which were more significant at low NaCl 
concentrations, resulted from the reduced dependence on the 
salt concentration, and this effect is more significant with large 
PEG molecules. 

2.3. Ribozyme reactions in other mixed solutions 

We also employed the following compounds that are inert or 
non-inert to RNA: ethylene glycol (EG, the monomer unit of 
PEG); glycerol (Glyc), 1,3-propanediol (PDO), 2-methoxyethanol 
(MME), and 1 ,2-dimethoxyethane (DME) (structurally related to 
EG): methanol (MeOH), ethanol (EtOH), and 1-propanol (PrOH) 
(small primary alcohols): urea and formamide (FA) (amide com- 
pounds, known as nucleic acid denaturants): N-methylformamide 
(NMF), dimethylformamide (DMF), and acetamide (AcAm) (struc- 
turally related to FA): acetonitrile (AcCN), dimethyl sulfoxide 
(DMSO), tetrahydrofuran (THE), and 1,4-dioxane (DOX) (aprotic 
solvent molecules): and polysaccharide dextran (Dex) with an 
average molecular weight 1 x 10*. Many of these compounds 
changed the cleavage rate (Fig. SIB), giving the rate constant range 
of 0.0029-3.7 h"' with 1.5 M NaCl, as compared in Fig. 2. In partic- 
ular, urea and FA significantly decreased the rate constant by 
100-fold, consistent with their role as nucleic acid denaturants. 
The FA-like compounds of NMF, DMF, and AcAm also decreased 
the value by 3.9- to 6.5-fold, whereas many of the other com- 
pounds increased or had slight effects on the rate constant. Some 
of the cosolute effects are consistent with prior studies measuring 
the hammerhead ribozyme activity with Mg^* in the mixed solu- 
tions with Glyc or MeOH [31 [. Dex as a polymer molecule had a 
lower influence on the cleavage rate, but the small primary 
alcohols that do not cause a large steric effect increased the rate 
constant as much as large PEG molecules did. The same effects of 
these compounds on the ribozyme catalysis were observed when 
the data at low NaCl concentrations were compared. The observa- 
tions suggest that factors besides the excluded volume effect have 
major effects. It is a remarkable finding that these cosolute mole- 
cules changed the slope of the linear plot of log k versus log [NaCl] 
(Fig. SIC), as summarized in Table 1. The slope values had a weak 
correlation with the k values at a certain NaCl concentration, e.g., 
1.5 M NaCl (Fig. SID). 
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Fig. 1. (A) RNA sequences of the hammerhead ribozyme and the FAM-labeled substrate RNA used in this study. (B) The polyacrylamide gel electrophoresis for the ribozyme 
reaction for 30 min with 10 mM MgCl2 or 0, 0.5, 1.0, or 1.5 M NaCI. The arrowheads indicate the uncleaved (upper) and cleaved (lower) fragments of the substrate RNA. The 
reactions with NaCI in the presence of PEG200 (P2h) or PEG8000 (P8k) at 20 wt% are also presented. (C) Changes in the reaction rate constant by the addition of PEG of 
different molecular weights, expressed as the ratio of the rate constants with (((pec) and without PEG (fcnone), in the solutions containing 1.5 M (circles) or 0.5 M NaCI 
(triangles). (D) The NaCI concentration dependence of k in the absence (black) and presence of PEG200 (blue), PEG600 (cyan), PEG2000 (orange), PEG8000 (red), or PEG20000 
(green) at 20 wt%. The correlation coefficients of a linear fit are greater than 0.995. (For interpretation of the references to color in this figure legend, the reader is referred to 
the web version of this article.) 
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Fig. 2. Rate constants for the ribozyme reaction in the 1 .5 M NaCI solutions without and with cosolutes at 20 wt%. 



Table 1 

Dependence of the rate constant of the ribozyme-catalyzed RNA cleavage on the NaCI 
concentration. 
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Values represent the slope of the plot of log k versus log [NaCI] obtained using the 
mixed solutions with a cosolute at 20 wt%. Each error value is within 0.2. 



2.4. Stability of RNA base pairs in the mixed solutions 

We also investigated tiie thermodynamic stability of RNA base 
pairs in the mixed aqueous solutions. In general, the base-pair for- 
mation is accompanied by the association of metal ions, and the 
ions preferentially interact with the phosphate groups through 



the Coulomb interaction, without obvious site specificity [32]. 
We prepared an ll-mer RNA duplex 5'-AUCCAGCGCGC-375'-GU 
GCGCUGGAU-3', derived from the base-pair sequence formed by 
the hammerhead ribozyme and the substrate RNA. Circular dichro- 
ism (CD) spectra of the duplex displayed no obvious change in the 
signal in the mixed solutions (Fig. S2A), suggesting that the double- 
helical structure is not significantly affected by the cosolutes 
although the spectra are not particularly sensitive to subtle or local 
conformational changes. 

The values of the free energy change for the duplex formation 
determined by UV melting curve analysis showed that no coso- 
lutes, excluding those with PEG8000 (phase separation at high 
temperatures), PEG20000 (phase separation at high temperatures), 
and THF (evaporated at high temperatures) for which measure- 
ments could not be obtained, increased the duplex stability at 
1 M NaCI (Fig. S2B). As found in the ribozyme study, the amide 
compounds had the strongest effect, e.g., FA decreased the Tm by 
9 °C and - AG° by 2.8 kcal mol^' (Fig. 3A), consistent with previous 
reports that the stability of DNA duplexes was decreased in the 
presence of FA or urea due to interactions with unpaired bases 
through the formation of multiple hydrogen bonds [33-35]. EG, 
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Fig. 3. (A) Melting curves of the 1 l-mer RNA duplex in the 1 M NaCI solutions without (black) and with PEG200 (blue), EtOH (brown), FA (red), or DMF (cyan) at 20 wt%. (B) 
Plots of log K versus log [NaCI]. The data for PEG8000 (purple) at high NaCI concentrations were not available due to phase separation. The correlation coefficients of a linear 
fit are greater than 0.990. (C) Plots between the NaCI concentration dependence of log k and that of log K. The data using PEG and EG are indicated by black symbols, and those 
using the amide compounds and other compounds are indicated in green and red, respectively. The data using the solution without cosolutes are indicated in gray. (For 
interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 



EG-like compounds (Glyc, PDO, MME, and DME), primary alcohols 
(MeOH, EtOH, and PrOH), and aprotic compounds (AcCN, DMSO, 
THE, and DOX), which increased or had slight effects on the RNA 
cleavage rate, significantly decreased the value of -AG°. Small 
PEG molecules decreased the duplex stability (1.6 kcal mol"' by 
PEG200 and 0.9 kcal mol^^ by PEG600), which can be explained 
by specific interactions with single-stranded bases [13] and reduc- 
tion in the water activity [14], Small cosolutes may act as osmo- 
lytes that reduce the activity of water and increase the osmotic 
pressure [36,37]. The solutions containing PEG or the amide com- 
pounds have moderate water activities (0.94-0.96), and those con- 
taining EG, the EG-like compounds, the small primaiy alcohols, and 
the aprotic compounds have relatively low values (0.93 or less). 
We found a weak correlation of the -AC°, excluding those of the 
amide compounds, with the logarithm of the water activity, while 
no correlation was observed with the viscosity (Fig. S3A and B). In 
addition, the k values of the ribozyme reaction showed no clear 
correlation with these solution property parameters (Fig. S3C and 
D), and there was no strong correlation between the values of k 
and -AC" (Fig. S3E). 

Some organic compounds may cause microphase separation or 
decrease the solubility of NaCI, which could increase the effective 
concentration of NaCI. The excluded volume effect of large PEG 
molecules also increases the effective concentration of NaCI, and 
the effect becomes more significant when a large PEG is used. These 
effects on the salt concentration are expected to enhance the duplex 
stability; however the studied mixed solutions did not increase the 
stability of the RNA duplex, excepting those containing large PEG 
molecules at low concentrations of NaCI. Fig. 3B shows the compar- 
ison of the effect of the duplex stability on the NaCI concentration. 
Linear dependences between log K (= - AG°/2.303RT, where R is the 
gas constant and lis the absolute temperature) and log [NaCI] with 
a correlation coefficient greater than 0.990 were obtained, where 
the data for the small primary alcohols at high salt concentrations 
and those for large PEG molecules causing phase separation at high 
salt concentrations and high temperatures were omitted from the 
linear regression analysis. The slope value 1.8 in the absence of 
cosolutes was within those previously obtained for DNA and RNA 
structures, e.g., ranging from 0.6 to 4 which varied depending on 
the length and structure [19,20,38]. The slope changed when using 
the mixed solutions: for example, the slope value decreased to 1.1 
in the PEG8000-containing solution. The reduction in the depen- 
dence on the salt concentration was also reported in previous stud- 
ies using DNA duplexes with PEG [ll,18-20[ and RNA tertiary 
structures with trimethylamine oxide (TMAO) [38]. Our results 
show that many of the mixed solutions changed the magnitude of 
the salt concentration dependence of the duplex stability. In 



contrast, the slope value did not decrease in the FA- or urea-con- 
taining solutions. This property is important to keep their role as 
nucleic acid denaturants at low salt concentrations. Taken together 
with the ribozyme data, there is found a correlation between the 
NaCI concentration dependence of log fc and that of log K (Fig. 3C) 
even though the ribozyme activity depends not only on the stability 
of secondary structures but also on tertiary interactions. 

2.5. Correlation with the dielectric constant 

We then considered the relative dielectric constant, Ep An 
inclusion of an organic compound having a low dielectric constant 
lowers the dielectric constant of a solution ] 39,40]. Most of the 
studied compounds decrease the dielectric constant, but some 
amide compounds have a unique dielectric property: the addition 
of FA or NMF increases the value of while DMF decreases the 
value. The studied solutions have values ranging from 55 to 100 
(Fig. S4). Although the creation of a uniform dielectric medium is 
assumed and the dielectric constant decreases with the salt 
concentration and the temperature [41 ], there is some correlation 
between the values of e^' and k of the ribozyme reaction (Fig. 4A). 
More strikingly, the inverse of e,- correlated well with the NaCI con- 
centration-dependent data given in Table 1: the slope value 
decreases as e^^ increases, except in the case of NMF 
(e^^ = 0.010), as shown in Fig. 4B. The deviation for the NMF-con- 
taining solution having the highest dielectric constant might be 
due to the overestimation of £r caused by a nonhomogeneous dis- 
tribution in the solution. Furthermore, the x-intercept of the log k 
versus log [NaCI] plots also showed a correlation with the e^' 
(Fig. 4C). It is emphasized that the data points obtained using the 
mixed solutions with amide compounds, but excluding NMF, fell 
along the same lines of the correlation plots. The data obtained 
using the solutions with different cosolute amounts ranging from 
5 to 20 wt% and the binary mixture solutions of PEG200, 
PEG8000, EtOH, and FA followed the same lines of the plots, but 
no correlation was observed with the water activity or viscosity 
(Fig. S5). The same correlations with the c^^ were also found for 
the duplex stability regardless of the cosolute property of interact- 
ing with RNA, but no correlation was observed between the values 
of -AG° and e^' (Fig. 4D-F). 

3. Discussion 

3.1. Influence of the mixed solutions on the ribozyme catalytic activity 

This study investigated the RNA cleavage rate of the hammer- 
head ribozyme and the thermodynamic stability of RNA base pairs 
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Fig. 4. (A) Plots of k obtained with 1.5 M NaCl versus of the solutions. Symbols are the same as in Fig. 3C. (B and C) Plots of the slope orx-intercept of the log k versus log 
[NaCl] plot against e^'. Symbols are the same as in Fig. 3C, and the data using the solutions with different cosolute amounts (5 or 10 wt%, given in Fig. S5) and the binary 
mixture solutions of PEG200 and PEG8000, PEG200 and EtOH, PEG200 and FA, or EtOH and FA, each at 1 0 wt%, are indicated in blue. The correlation coefficients of a linear fit, 
excluding the data using NMF (e^' = 0.01 0), are 0.861 in panel B and 0.772 in panel C. (D) Plots between -AC of the RNA duplex at 1 M NaCl and e^'. The data using PEG8000, 
PEG20000, and THE are excluded because of the phase separation or evaporation at high temperature during thermal melting. Symbols are the same as in Fig. 3C. (E and F) 
Plots of the slope orx-intercept of the log /C versus log [NaCl] plot against e,^'. The correlation coefficients of a linear fit, excluding the data using NMF, are 0.906 in panel E and 
0.888 in panel F. Symbols are the same as in panels B and C. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this 
article.) 



in the solutions containing NaCl. The experiments were performed 
with Na* rather than Mg^* because of several benefits. First, 
monovalent ions are less likely to associate with cosolute mole- 
cules. Second, the association of Na* ions is mostly nonspecific 
through electrostatic interactions. The ion binding enhances the 
activity of the hammerhead ribozyme by increasing the population 
of the catalytically active conformation relative to inactive confor- 
mations. Third, Na* ions weakly but cooperatively associate with 
RNA, which can be analyzed by a linear regression method pre- 
sented in Figs. ID and 3B. The plot of the NaCl concentration 
dependence provides a steeper slope than the case of MgCb 
because a greater number of monovalent ions are required for 
the charge neutralization than when using divalent ions [42]. This 
facilitates the comparison of the salt concentration-dependent 
data obtained using different solutions. 

The studies on the effect of PEG on the hammerhead ribozyme 
showed increased rates of RNA cleavage in the solutions containing 
large PEG molecules. The results are consistent with the excluded 
volume effect that may increase the effective concentration and 
restrict the degree of freedom of RNA molecules [10]. However, 
the k value did not continuously increase with the molecular 
weight of PEG. The result suggests that the simple excluded vol- 
ume model cannot fully explain the PEG effects. The k values also 
have no clear correlation with the water activity, viscosity, or 
RNA base-pair stability of the ll-mer duplex, but have some corre- 
lation with the r,^^ (Fig. 4A). A simple interpretation of the correla- 
tion between the rate constant and the dielectric constant is that 
the transition state is electrostatically stabilized in less-polar 
media. However, since the rate of RNA cleavage with monovalent 
ions is controlled by the folding process rather than by the catalytic 
step, it is more likely that the low dielectric environment enhances 
the folding into the catalytically active conformation mediated by 
sodium ion binding (discussed later). It has been reported that 



the RNA cleavage rates of the hairpin ribozyme and the lead- 
dependent ribozyme (leadzyme) are accelerated in the mixed solu- 
tions with PEG (MW = 4x 10^): the acceleration of the hairpin 
ribozyme reaction with Mg^* was attributed to the stabilization 
of the catalytically active conformation under the PEG-mediated 
osmotic pressure, and not to the influence of the dielectric prop- 
erty ]43,44]. The acceleration of the leadzyme reaction with Pb^* 
was also attributed to the osmotic pressure effects by PEG [45]. 
On the other hand, the study for the group 1 intron ribozyme motif 
showed that the PEG (MW = 1 x 10^)-containing solution favored a 
more compact structure of the large RNA molecule and enabled the 
folding with lower Mg^* concentrations, attributed to the domi- 
nant influence of the excluded volume effects [46[. It is possible 
that the origin of the PEG effect varies among the different ribo- 
zyme motifs, due to the differences in the structure and stability 
of RNA, molecular sizes of RNA and PEG, rate-limiting step, and 
contribution of ions and waters to the reactions. 

3.2. Prediction method in the mixed solutions 

It is convenient for practical use to predict nucleic acid struc- 
tures and their stability under non-dilute conditions. The thermo- 
dynamic parameters based on the nearest-neighbor interaction 
model can be used for the predictions of the hybridization energy 
and secondary structures of nucleic acids [15[, and we have previ- 
ously found the applicability of the nearest-neighbor model for 
DNA structure formations in several mixed solutions [14[. It was 
reported that EG, MeOH, EtOH, PrOH, urea, FA, DMF, DMSG, and 
DOX destabilized a short RNA duplex, and there was found no good 
correlation between the duplex stability and the viscosity, dipole 
moment, surface tension, or solubility of the nucleotide adenine 
[47]. We demonstrated in this study that there was no correlation 
between the -AG° of the ll-mer RNA duplex and the dielectric 



648 



S. Nakano et al./FEBS Open Bio 4 (2014) 643-650 



constant (Fig. 4D). These results are probably due to significant 
energetic contributions from specific interactions witli several 
cosolutes {e.g., amide compounds) and from base pairing in a rela- 
tively low water activity environment {e.g., EG, Glyc, and small pri- 
mary alcohols). 

The dependence of the duplex stability on the NaCl concentra- 
tion was characterized by a linear correlation between log K and 
log [NaCl], as shown in Fig. 3B. Although increasing the ionic 
strength may replace cosolute molecules that possibly accumulate 
at the nucleic acid surface, the linear correlations suggest that the 
replacements could have only minor effects on the data analysis. 
To evaluate the salt concentration dependence, empirical expres- 
sions that take into account the number of phosphate charges 
contributing to base pairing and the Coulombic end effects in ion 
accumulation have been developed [16,17]. Our results demon- 
strate that the salt concentration dependence changes as a function 
of the dielectric constant (Fig. 4E), which could have some rele- 
vance to the phosphate charge effect and the Coulombic end effect. 
Based on the finding, the value of -AG° of an RNA duplex under 
non-dilute conditions may be predicted by considering a possibil- 
ity of specific interactions with coexisting compounds, water activ- 
ity effects, and dielectric constant effects, for example, by using the 
following calculation: -AG° at a certain salt concentration in a 
non-dilute solution = (-AG° at a standard condition, e.g., 1 M NaCl 
in a dilute solution) + (energy penalty due to specific interactions 
and the reduced water activity) + (a term involving the dielectric 
constant, possibly combined with the phosphate charge effect 
and the end effect) x log[salt]. More comprehensive studies using 
other sequences and concerning the energy penalty and sequence- 
and length-dependent dielectric constant effects will contribute to 
nucleic acid folding predictions under non-dilute conditions. 

3.3. Dielectric constant effects on nucleic acids 

Based on the polyelectrolyte model of nucleic acids, the magni- 
tude of the dependence of duplex stability on the salt concentra- 
tion is related to the thermodynamic degree of ion accumulation 
near nucleic acids 148-50]. Our finding that mixed solutions chan- 
ged the NaCl concentration dependence (the slope of the linear 
plots in Fig. 3B) suggests altered binding properties of salt ions to 
RNA molecules in the solution conditions created by organic com- 
pounds. It can therefore be assumed that the linear correlations 
between the salt concentration dependence data and the c^^ 
indicate the significance of electrostatic ion binding to RNA, med- 
iated by Coulomb interaction that is inversely proportional to the 
dielectric constant. In addition to the slope of the linear plots, their 
x-intercept values have a correlation with the c^' (Fig. 4C and F). 
The environment with a low dielectric constant would enhance 
the efficiency of ion binding to RNA molecules, and is favorable 
both for the ribozyme catalysis and for the duplex formation with 
a lower amount of NaCl, and one of such compounds is PEG having 
a molecular weight of more than several thousands. In contrast, 
urea and FA providing a medium of high dielectric constant effec- 
tively destabilize the RNA structures in both high and low salt con- 
centrations. Their role as nucleic acid denaturants could be 
attributed not only to the formation of hydrogen bonding with sin- 
gle-stranded bases, but also to the property of not creating a low 
dielectric constant condition. 

It is noted that the values of the dielectric constant may not 
accurately reflect the local environment near the RNA surface. 
However, it is reported that ions distant from an RNA surface inter- 
act significantly with an RNA, and long-range interactions are 
suggested to be important in ion-induced stabilization for RNA 
structures ]51 ]. Our observations imply that ions at some distance 
from the RNA surface, influenced by the bulk dielectric constant, 
play an important role in the stabilization of the RNA structures: 



sodium ions are distributed around the phosphate backbone and 
strong electronegative positions in the active site for both the reac- 
tant and activated precursor states of the hammerhead ribozyme 
127], and a number of sodium ions bind primarily to the phosphate 
groups of a base-paired duplex 132]. Importance of the dielectric 
constant on nucleic acid reactions involving metal ion binding 
has also been proposed for the multivalent ion-mediated DNA 
condensation and aggregation 152,53], thermal stability of a DNA 
G-quadruplex structure formed by coordination of K* ions 154], 
structural changes in a functionally important bulge region of the 
group II intron ribozyme and its Mg^* binding properties 155], 
and dimerization of RNA or DNA hairpins induced by metal ion 
binding [56]. A relatively low dielectric constant of about 50 or 
even less, which is significantly lower than the value for pure water 
of approximately 80, is suggested for the inside of intact cells 
157,58]. A computer simulation study predicts a reduced dielectric 
constant values in the range of 30-50 under the molecular crowd- 
ing condition 159]. In addition to the intracellular environment, the 
mobility and arrangements of water near the surface of nucleic 
acid biosensors are highly restricted, and thus the dielectric con- 
stant of the medium of a nucleic acid probe is significantly lower 
than that of a dilute solution. There would be important contribu- 
tions of the low dielectric constant environment inside of cells and 
near the surface of a biosensor to the ribozyme activity and the sta- 
bility of nucleic acid structures. The molecular crowding effects on 
nucleic acids are well explained by the excluded volume effect, the 
osmotic pressure effect, or specific interactions, and it would also 
be important to consider the effects of the change in the dielectric 
constant in the molecular crowding studies using mixed solutions 
with organic compounds. 

4. Conclusions 

We have found that many types of mixed solutions change the 
rate of the ribozyme-catalyzed RNA cleavage and its dependence 
on the concentration of NaCl. The solutions also affect the stability 
of RNA base pairs and its dependence on the salt concentration. 
Specific interactions with RNA molecules (e.g., urea and FA) and 
the reduced water activity (e.g., EG, Glyc, and small primary alco- 
hols) can explain their inhibiting effects on the ribozyme catalysis 
and decreasing effects on the duplex stability, but they do not 
explain the changes in their NaCl concentration dependence. We 
found good correlations of the dielectric constant with the salt con- 
centration dependence of the ribozyme reaction rate and the RNA 
base-pair stability. The observation provides useful insights when 
constructing a prediction model for RNA structures under 
non-dilute conditions and also when modulating the stability by 
changing the solution composition. 

5. Materials and methods 

5.1. Preparation of RNA oligonucleotides and buffer solutions 

The hammerhead ribozyme sequence derived from Schistosoma 
mansoni was enzymatically prepared as previously described 130]. 
The substrate RNA labeled with 6-carboxyfluorescein (FAM) at the 
5'-end and the 11-mer RNA oligonucleotides, purified by high-per- 
formance liquid chromatography (HPLC), were purchased from 
Hokkaido System Science. 

All reagents used to prepare buffer solutions were purchased 
from Wako with the following exceptions: 4-(2-hydroxyethyl)-l- 
piperazineethanesulfonic acid (HEPES) and the disodium salt of eth- 
ylenediamine-N,N,N',N'-tetraacetic acid (Na2EDTA) from Dojindo, 
PEG8000 and dextran with an average molecular weight 1 x 10"* 
from Sigma, and 2-methoxyethanol and 1,2-dimethoxyethane from 



S. Nakano et at./FEBS Open Bio 4 (2014) 643-650 



649 



TCI. All reagents were used without further purification. Cosolute 
molecules were added to buffer solutions at 20 wt%, unless 
otherwise mentioned. No obvious precipitation and phase separa- 
tion were observed, excluding the experiments using large PEG 
molecules at high salt concentrations and high temperatures. 

5.2. Ribozyme-catatyzed RNA cleavage 

The RNA cleavage rate of the hammerhead ribozyme was mea- 
sured under the single-turnover conditions, in which the ribozyme 
(2 iiM) was present in excess relative to the substrate (100 nM). 
The reaction was performed at 37 °C using a buffer solution 
consisting of 50 mM HEPES and 1 mM Na2EDTA at pH 7.0. No sig- 
nificant contamination of divalent metal ions was confirmed from 
the experiments using higher concentrations of Na2EDTA. 

The RNA solution was annealed in the buffer solution at 60 °C 
and then incubated at 37 °C for 10 min, followed by the addition 
of salts to start the reaction, during which the addition of NaCl 
changed the solution pH by no more than 0.13. The quenched 
reaction mixture was loaded onto a 20% polyacrylamide gel (acryl- 
amideibisaciylamide = 19:1) containing 7M urea. After the gel 
electrophoresis, the fluorescence emission from the FAM-labeled 
RNA strands was visualized and quantified by a fluorescent scanner 
(FLA-5100, Fuji) using a 473-nm excitation laser and a 520-nm 
emission filter. 

The observed rate constant, k, for the ribozyme-catalyzed cleav- 
age of the substrate RNA was determined from a plot of the RNA 
fraction cleaved at the correct site,/, versus time, t, fitted to the sin- 
gle-exponential equation of/=/o + (/max -/o) [1 - exp(-/<:t)], where 
/o is the fraction at time zero and/^ax is the fraction at the endpoint 
of the reaction. The fitting was performed using KaleidaGraph 4.1 
(Synergy). A linear approximation of the equation was used to 
calculate the rate constant of very slow reactions. 

5.3. CD spectra of an RNA duplex 

CD spectra were obtained using a spectropolarimeter (J-820, 
JASCO). All spectra were measured at 5 °C using an RNA strand con- 
centration of 40 |.iM in a buffer consisting of 1 M NaCl, 10 mM Na2- 
HPO4, and 1 mM NaiEDTA at pH 7.0. The RNA solution was heated 
to 80 °C and cooled at the rate of 2 °C min"' prior to use. 

5.4. Thermodynamic stability of an RNA duplex 

The thermal melting curve of the RNA duplex was obtained by 
monitoring the absorption at 260 nm using a spectrophotometer 
(UV-1800, Shimadzu). The melting curve was measured using an 
RNA strand concentration of 2 \iM in the phosphate buffer. The 
heating rate was 1 °C min"\ and the cuvette was sealed with an 
adhesive sheet to prevent evaporation. The melting temperature. 
In,, at which half of the duplex structure was denatured, the asso- 
ciation equilibrium constant, K, and the free energy change, AG°, at 
37 °C were determined from the melting curve [60]. 

5.5. Properties of the mixed aqueous solutions 

The water activity was determined by the osmotic pressure 
method using vapor phase osmometry (5520XR pressure osmom- 
eter, Wescor) or freezing point depression osmometry (Typ Dig.L 
osmometer, Halbmikro). The relative dielectric constant was 
calculated according to the equation reported by Oster or experi- 
mentally determined using a fluorescent probe, l-anilino-8-naph- 
thalene sulfonate, as previously described [56]. The solution 
viscosity was measured with a viscometer (SV-10 vibro viscome- 
ter, A&D). These solution property parameters were measured 



using the buffer solutions containing 1 M NaCl and 20 wt% cosolute 
at 25 °C, unless otherwise mentioned. 
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